Heavy Quark Production from Jet Conversions in a Quark-Gluon Plasma 
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Recently, it has been demonstrated that the chemical composition of jets in heavy ion collisions is 
significantly altered compared to jets in the vacuum. This signal can be used to probe the medium 
formed in nuclear collisions. In this study we investigate the possibility that fast light quarks and 
gluons can convert to heavy quarks when passing through a quark gluon plasma. We study the 
rate of light to heavy jet conversions in a consistent Fokker-Planck framework and investigate their 
impact on the production of high-pr charm and bottom quarks at RHIC and LHC. 
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I. INTRODUCTION 

In heavy ion collisions, energetic jets are expected to 
lose a significant amount of energy when passing through 
the surrounding quark gluon plasma (QGP) [J [3, H, [j| • 
This jet quenching phenomenon leads to a large suppres- 
sion of the yield of high-py hadrons in Au+Au collisions 
compared to theyield expected from a superposition of 
p + p collisions 0, la] at the same center-of-mass energy. 
Induced gluon radiation was expected to be the largest 
contributor to the momentum degradation of the lead- 
ing jet parton. However, experimental results on the 
quenching of heavy flavors Ml has led to new efforts to 
revisit elastic energy loss [1,1a, G3 and non-perturbative 
mechanisms [Tlj . Jet quenching measurements deliver 
important information about the medium in terms of the 
average squared momentum transfer to the jet per unit 
path length, q = /i 2 /A. Besides quenching as a manifes- 
tation of energy loss of the leading high-pr parton, the 
chemical composition of a jet can also change through fla- 
vor changes of the leading parton. This is an unavoidable 
consequence of collisions of the high-pr parton with ther- 
mal partons. The probability of conversion is related to 
the path length through the medium and the collisional 
conversion widths. Thus by tracking the flavor of jets 
(defined as the flavor of the leading parton) throughout 
their interaction with the medium, transport properties 
of the medium which are independent of and complemen- 
tary to q can be obtained. Eventually, the flavor changes 
will be reflected by the hadron composition of the jet. 

Flavor changing processes were first studied in the con- 
text of light quarks and gluons [l2j , and parton to pho- 
ton conversions jTJ . More recently, we unified these two 
processes and discussed new observables which could be 
measured at the Relativistic Heavy Ion Collider (RHIC) 
or the Large Hadron Collider (LHC) [H|. In Ref. [13 
it was found that conversions of light quarks to gluons 
could help solve the puzzle of very similar nuclear mod- 
ification factors Raa for pions and protons observed by 
the STAR experiment [14(. These seemed to be incom- 
patible with a relative energy loss of 9/4 for gluons and 
light quarks. However, as we discussed in more detail in 
[La ] , the concept of a fixed flavor for the leading jet parton 



is ill-defined in a medium. Conversions to real or virtual 
photons are just other examples of what can happen to 
a fast quark or gluon traversing a quark gluon plasma. 
Although much rarer than conversions to strongly inter- 
acting partons, it has been shown that these photons can 
make a significant contribution to the direct photon or 
dilepton spectrum at intermediate px- 

In [lH it was shown that strange hadrons might be the 
ideal probe for flavor conversions at RHIC energies. This 
is due to the strong gradient between the very small ini- 
tial strange quark content of jets at RHIC, and the almost 
full chemical equilibration of strangeness in the plasma. 
Jets coupling to the medium will be driven toward chemi- 
cal equilibrium and thus kaons and A-hyperons have to be 
enhanced relative to pions and protons at large px when 
compared to naive expectations from p + p collisions. 
Measurements of this enhancement could shed light on 
the mean free path of fast gluons and light quarks in the 
medium. At LHC, the value of strangeness as a high-pr 
probe is diminished through the large initial strangeness 
content of jets at larger energies. When translating lead- 
ing particle flavor into hadrons it has to be kept in mind 
that the hadron chemistry can also be altered through 
the increased multiplicity in a jet cone interacting with 
the medium as studied in (l6| . 

Two important questions emerge from the discussion 
of high-pr strangeness enhancement at RHIC. First, can 
heavy quarks, in particular charm, play a similar role as 
a probe for the average mean free path at LHC? And sec- 
ondly, is the suppression of charm found at RHIC com- 
patible with the mechanism proposed here which boosts 
the yields of rare particle species at intermediate and 
high pt? In this Report, we answer these questions by 
studying the production of high-p^ heavy quarks from 
jet conversions in heavy ion collisions. We compute the 
contributions from leading order (LO) processes in per- 
turbation theory multiplied by a K- factor. The heavy 
quarks can be produced from the annihilation process 
g + g — ► Q + Q that converts a gluon jet and a thermal 
gluon to a pair of heavy quarks, and also through the 
Compton processes g(q) + Q — > Q + g(q) by transferring 
the momentum of the incoming jet to that of a slow mov- 
ing heavy quark in the medium. We follow the method 
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introduced in Ref. [15| to study the propagation of the 
high pt heavy quark in the expanding fireball within the 
framework of a Fokker-Planck equation. 

This work is organized as follows. In the next section 
we discuss the transport coefficients of heavy quarks and 
the collisional widths for the conversion of light partons 
to heavy quarks. We then proceed to present numerical 
results for charm and bottom quarks at RHIC and LHC 
in Sec lIIII using a consistent model based on rate equa- 
tions for conversions and a Fokker-Planck equation for 
energy loss. Finally, a summary is presented in Sec. [TV] 



II. HEAVY QUARK SCATTERING IN A QGP 

First we discuss the propagation of heavy quarks in a 
partonic medium. It can be described in the framework 
of a Fokker-Planck equation [ll|, [I?} 



df( P ,t) = d_ 
dt <9pi 
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f(p,t) (1) 



where /(p, t) is the distribution function of heavy quarks. 
The drag and diffusion coefficients are 
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Here p and p' label the momenta of the heavy quark 
before and after each collision, respectively. We define 
the thermal averaging (Y"(p')) in the elastic process Q + 
g(q) -> Q + g{q) as 
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where M. is the scattering amplitude and jq = 6 is the 
spin and color degeneracy factor for heavy quark. 

Using a fixed coupling constant a s = 0.3 and masses 
m g = <?T/\/2 for thermal gluons and m q = gTj\f% 
for thermal quarks jl8j ]. we calculate the transport co- 
efficients for both charm (m c =1.2 GeV) and bottom 
(rnft=4.5 GeV) quarks numerically at leading order ac- 
curacy for the Compton processes Q + g(q) — > Q + g(q). 
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FIG. 1: (Color online) Transport coefficients 7 of charm (left) 
and bottom quarks (right) in a QGP as functions of transverse 
momentum px at different temperatures T. 
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FIG. 2: (Color online) Transport coefficients Bo of charm 
(left) and bottom (right) quarks in a QGP as functions of 
transverse momentum pr at different temperatures T. 



We show the results in Figs. [TJ and [3] The drag co- 
efficients 7 decrease with increasing momentum of the 
incoming jets while diffusion coefficients behave the op- 
posite way. 

To study the dynamics of heavy quark in an expanding 
QGP medium, we focus on the degradation of transverse 
momentum of heavy quark jets, which can be taken into 
account using a Langevin equation (l9l [20| obtained from 
the Fokker-Planck equation in (fT]) by calculating the first 
and second modes, 
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Here AW is the random force according to the normal 
distribution 
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and Bij is the inverse of the diffusion-coefficient matrix 
Bij defined in Eq. ©. 
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FIG. 3: (Color online) Transport coefficients B\ of charm 
(left) and bottom (right) quarks in a QGP as functions of 
transverse momentum pr at different temperatures T. 



Now we discuss the implementation of the conver- 
sion of light quarks (u,d,s) and gluons into heavy quarks 
through the processes g + g — > Q + Q and g{q) + Q — > 
Q + g(q). The dynamics is governed by the rate equation 
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quark after conversion from a gluon, pq = (pq)/(1), as a 
function of the incoming gluon momentum p are shown 
in Fig. [5l {pq) has an almost linear dependence on p only 
depends very weakly on temperature. 
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FIG. 5: (Color online) Average final state momentum of the 
heavy quark for charm (left panel) and bottom (right panel) 
after conversion from gluons with momentum p for two dif- 
ferent temperatures. 



The collisional conversion widths can be estimated as 
T(p) = (1) in the notation denned in Eq. Q. 
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FIG. 4: (Color online) Conversion widths for charm (left 
panel) and bottom (right panel) in a QGP due to the pro- 
cess g + g^Q + Qas functions of transverse momentum pr 
of the incoming gluon jets at different temperatures T. 

In Fig. IU we plot the width of a gluon jet converting 
to charm and bottom via the process g + g — > Q + Q. 
We have introduced a constraint on the momentum of 
the outgoing heavy quark jet. We only count the heavy 
quark as a jet particle if it carries the larger momentum 
of the two final state partons. The conversion will also 
cause effective energy loss since the final momentum is 
usually smaller than the initial gluon jet momentum due 
to the mass threshold. This is different from the case 
of light flavor conversions which can occur with 100% 
efficiency. The resulting average momentum of the heavy 



As mentioned above there are two more elastic pro- 
cesses, g+Q — > Q+g and q+Q — > q+Q, which contribute 
to the enhancement of high-pT heavy quarks by acceler- 
ating slowly moving heavy quarks from the medium. We 
use the distribution function of heavy quarks obtained 
from perturbative QCD calculations in Pb+Pb collisions 
at TeV (shown in the next section) to esti- 

mate the conversion width for a quark or gluon into a 
heavy quark via these Compton-like processes. It turns 
out that their contribution to high-py heavy quark pro- 
duction is two orders of magnitude smaller than that from 
the gluon annihilation process g + g — ► Q + Q due to the 
mass mismatch in the initial state. This is true when we 
require the momentum of the outgoing heavy quark to be 
larger than half of the momentum of the incoming light 
flavor jet. When this constraint is lifted a much larger 
conversion rate is obtained, however, with a much smaller 
average momentum gain for the final state heavy quark 
which seems irrelevant for the high-py sector. Thus we 
neglect the contributions from Compton-like elastic pro- 
cesses in the following. 

In this work we follow the procedure outlined for light 
quark and gluons in [121 ]. The leading order matrix el- 
ements are scaled by a factor K = 4 which was needed 
to reproduce the observed p/n + and p/n" ratios. Ulti- 
mately, however, the drag and diffusion coefficients have 
to be extracted from data and compared to a variety of 
perturbative and non-perturbative predictions. A large 
iT-factor as used here, if found compatible with data, 
might rather hint toward a much stronger jet-medium 
coupling than expected from perturbation theory. 



4 



III. HEAVY QUARK CONVERSIONS AT RHIC 
AND LHC 

Now we turn to charm and bottom production in 
Au+Au collisions at RHIC with a center-of-mass energy 
yj snn — 200 GeV. The initial pt spectra of charm and 
bottom quarks at mid-rapidity are taken from Ref. [2l| . 
Both are obtained by multiplying the heavy quark px- 
spectra from p + p collisions at the same energy by the 
number of binary collisions (~ 960) in central Au+Au 
collisions. 
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FIG. 6: (Color online) The ratio of heavy quarks from jet con- 
versions to heavy quarks from initial production with energy 
loss included for charm (left panel) and bottom (right panel) 
as functions of transverse momentum pr in Au+Au collisions 
at ^s7«v = 200 GeV. 

For the dynamics of the fireball, we assume that it 
evolves boost invariantly in the longitudinal direction 
but with an accelerated transverse expansion. Specifi- 
cally, its volume expands in the proper time r according 
to V{t) = ttR(t) 2 t, where R(t) = R + o/2(t - r ) 2 
is the transverse radius with an initial value Rq=7 fm, 
To=0.6 fm is the QGP formation time, and a — 0.1c 2 /fm 
is the transverse acceleration [l2|, HJ, [22[ . With an initial 
temperature T, ; = 350 MeV and using thermal masses for 
quarks and gluons, this model gives a total transverse en- 
ergy comparable to that measured in experiments. The 
time dependence of the temperature is then obtained 
from entropy conservation, and the critical temperature 
T c = 175 MeV is reached at proper time r c ~ 5 fm. 

To solve Eq. [5J we follow the test particle methods in- 
troduced in Ref. for the propagation of light flavored 
jets and heavy quarks in the QGP medium. The light 
flavors can convert into each other or into heavy quarks 
according to the probabilities discussed above in [15j |. 
The conversion of heavy quarks into light flavors is sup- 
pressed and not taken into account, except for the process 
Q9(o) ~ * 9(q)Q which is included in the drag coefficient. 
In Fig. [5] we show the ratio of heavy quarks obtained 
from jet conversions to that obtained from initial hard 
collisions with energy loss included. The jet conversions 
contribute about 10% to the spectrum of charm quarks 
at high px, and only 0.3% to that of bottom quarks in 



central Au+Au collisions at RHIC. 

For LHC energies we take the initial pr spectra of 
charm and bottom quarks at midrapidity from the per- 
turbative calculation in Refs. [HI, [24| multiplied by the 
number of binary collisions ((N co u) sa 1700 [25]). The 
spectra are parametrized as 
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The parameters of the fireball are taken from Ref. [26j 
where thermal charm production at LHC was studied. 
Specifically, we choose initial proper time tq — 0.2 fm/c 
for the formation of the equilibrated QGP with initial 
temperature To = 700 MeV. In Fig. we show the ratio 
of heavy quark yields from jet conversions to that from 
initial hard collisions with energy loss included. The ra- 
tios are similar in magnitude magnitude to what has been 
obtained for RHIC. 



a 



o 



0.12 
0.09 

. 0.06 k, 

i 
) 

0.03 



— dN c /dN c 

COIW I 

K = 4 



Pb+Pb @ 5.5 TeV 



5 10 15 20 25 30 35 40 

p T (GeV/c) 



0.010 
0.008 

f 0.006 

v > 

8 0.004 

3 

0.002 
0. 



— dN b /dN b 

conv iro 

K = 4 



5 10 15 20 25 30 35 40 

p T (GeV/c) 



FIG. 7: (Color online) The ratio of heavy quark spectrum 
from initial production to that from jet conversions for charm 
(left panel) and bottom (right panel) as functions of trans- 
verse momentum in Pb+Pb collisions at^/Sjvjv = 5.5 TeV. 



We can now answer the two questions that we have 
posed at the beginning. First, the promotion of heavy 
quarks from the medium to high p? is not an effective 
mechanism for the production of high-p^ charm and bot- 
tom quarks. This is compatible with the experimental 
observations Q • Note that we have used the same frame- 
work and A"-factors which led to the prediction of a siz- 
able strange quark excess at high-pr at RHIC. On the 
other hand, we have found that the corresponding pro- 
cess for charm quarks at LHC is not efficient. This can 
be traced back to the relatively large amount of initial 
charm in jets produced at LHC and the missing chemical 
equilibration of charm in the QGP. Note that some au- 
thors have argued recently that thermal charm produc- 
tion might occur at the temperatures reached at LHC 
1201. We have chosen not to include this scenario here. 



5 



IV. DISCUSSION 

We have studied the impact of jet conversions on the 
production of high-p^ heavy quarks in heavy ion colli- 
sions. Based on leading order perturbative QCD, we cal- 
culated transport coefficients of heavy quarks, and stud- 
ied their propagation in an expanding QGP medium. We 
find that the dominant conversion process is pair pro- 
duction off gluon jets in the medium, g + g — ► Q + Q. 
However, the contribution of jet conversions to the total 
yields of heavy quarks at high px is rather small both at 
RHIC and LHC. This is consistent with measurements of 



single electrons at moderate values of pr at RHIC. On 
the other hand, it implies that high-py charm quarks at 
LHC as a probe for the mean free path in the medium 
might not be as useful as strangeness is at RHIC energies. 
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